**FULL TITLE** 

ASP Conference Series, Vol. **VOLUME** **YEAR OF PUBLICATION** 
**NAMES OF EDITORS** 



3D photoionisation and dust RT modelling with MOCASSIN: 
geometry effects on the emission line spectra of star-forming regions 



Barbara Ercolano & Nate Bastian 

Dept of Physics and Astronomy, University College London, UK 

Abstract. Emission line spectra from H II regions are often used to study prop- 
erties of the gas in star-forming regions, as well as temperatures and luminosities of 
the ionising sources. Empirical diagnostics for the interpretation of observational data 
must often be calibrated with the aid of photoionisation models. Most studies so far 
have been carried out by assuming spherical or plane-parallel geometries, with major 
limitations on allowed gas and dust density distributions and with the spatial distribu- 
tion of multiple, non-centrally-located ionising sources not being accounted for We 
present the first results of our theoretical study of geometric effects, via the construc- 
tion of a number of 3D photoionisation models using the MOCASSIN code for a variety 
of spatial configurations and ionisation sources. We compare integrated emission line 
spectra from such configurations and show evidence of systematic errors caused by the 
simplifying assumption of a single, central location for all ionising sources. 



Effects of 3D geometry on the emission line spectra of ionised regions 

The MOCASSIN code: Spatially-resolved studies of star-forming regions indicate that 
the assumption of spherical geometry is not realistic in most cases, further complication 
is added by the gas being ionised by multiple non-centrally located sources. The 3D 
photoionisation code MOCASSIN was designed to treat all of the above, allowing for the 
transfer of both primary and secondary radiation to be treated self-cons istently, without 
the need of approximations. The code was thoroughly benchmarked fercolano et al. I 
12003. 2005) and has been applied to the study of several ionised regions. The current 
version includes a fully self-consistent treatment of the radiative transfer of dust grains 
mixed within the gas, taking into account the microphysics of dust-gas interactions 
within the geometry-independent Monte Carlo transfer. 



3D gas and star distribution on the emission line spectra of ionised regions: We 

are investigating geometric effects theoretically, via the construction of a number of 3D 
photoionisation models for a variety of spatial configurations and ionisation sources. 
We compare integrated emission line spectra from such configurations in search of sys- 
tematic errors, which may be caused by the simplifying assumption of a single central 
location for all ionising sources. Our study is initially tailored at H II region spectra, 
but will successively be extended to star clusters and cluster complexes. 

Our preliminary models are based on homogeneous or porous spherical density 
distributions ionised by 240 sources, comprising of 37 Mq (35 kK) stars and 56 Mq 
(56 kK) stais. The stars are either located at the centre or distributed in the half- 
or full-spherical volum e. The ionising fluxes were calculated using STARBURST99 
jLeitherer et al. Ill999h . assuming solar metallicities and an age of IMyr. 
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Figure 1. Homogeneous sphere model; small, medium and large circles (P=0.72, 
0.74, 0.62) for central, half- and full-volume distributed sources. Left: Model R23 
ratios (circles) compared to Pilvugi^ ll2001h calibrations for P=0.2 (solid line), 0.35, 
0.5, 0.65 and 0.80. Right: Model ionic temperature (in [kK]) compared to assumed 
relations; solid line is for equal temperatures, dotted and dashed lines are the t[0 III], 
t[0 II] and t[0 III], t[S III] relations of lGarnettl(ll992l) . 



Preliminary results: R23 O-abundance diagnostic The sum of the bright O for- 
bidden lines ([O III] AA5007 ,4959 -i- [O ii]AA3727,29)/H/3 is often the sole handle on 
metallicity for extragalactic H II regions and starburst galaxies. We tested the robust- 
ness of abundances derived via this method using the calibration of Pilvugin (2001) . 
O abundance in our models for ce ntral half-volume distributed sources agree extremely 
well to the values obtained using "Pilvugin ' ("2001") calibrations for the appropriate ex- 
citation index, P. Significant deviations are shown, however, by the models with full- 
volume distributed sources For an homogeneous sphere (nH=100cm^'^), this implies 
that O abundances could be underestimated by ^25% by empirical methods (Fig 1, 
left). Similar results were found for the other density distributions investigated so far. 

Preliminary results: Ionic temperature relations We tested the robustness of rela- 
tions commonly assumed to hold between different ionic temperatures. Good agree- 
ment is found for central and half-volume stellar configuration cases in homogeneous 
density spheres, however models where the stars are distributed over the full volume 
have significantly lower temperatures than would be estimated from ID models (Fig. 1, 
left) for lower excitation species, such as 0+ and N+ and higher for higher excita- 
tion ones. This means that the t[0 III], t[0 ll] relation of Garnett ( 1992) fails for the 
full-volume stellar distribution models. Similar results were found for porous density 
distributions; here, however, larger variations are found for the t[0 III], t[0 ll] relation, 
with t[0 ll] severely underestimated for all stellar configurations. 
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